Random walks on fractals and stretched exponential relaxation 
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Stretched exponential relaxation (exp— (t/r) A ) is ob- 
served in a large variety of systems but has not been ex- 
plained so far. Studying random walks on percolation clus- 
ters in curved spaces whose dimensions range from 2 to 7, we 
show that the relaxation is accurately a stretched exponen- 
tial and is directly connected to the fractal nature of these 
clusters. Thus we find that in each dimension the decay ex- 
ponent fix is related to well-known exponents of the percola- 
tion theory in the corresponding flat space. We suggest that 
the stretched exponential behavior observed in many complex 
systems (polymers, colloids, glasses...) is due to the fractal 
character of their configuration space. 

PACS numbers: 61.20.Lc, 64.60.Ak, 05.20.-y 



The stretched exponential decay function, q(t) = 
exp — (t/r)^ K , was initially proposed empirically in 1854 
by R. Kohlrausch to parametrise the discharge of Ley- 
den jars, and was rediscovered in 1970 by Williams and 
Watts |Q . Since then the phenomenological " KWW" ex- 
pression has been shown to give an excellent represen- 
tation of experimental and numerical relaxation data in 
a huge variety of complex systems including polymers, 
colloids, glasses, spin glasses, and many more. This be- 
havior has attracted considerable curiosity; it has been 
discussed principally in terms of models where individual 
elements relax independently with an appropriate wide 
distribution of relaxation times (see for instance 
However, despite the ubiquity of the KWW expression, 
in the view of many scientists its status remains that of 
a convenient but mysterious phenomenological approxi- 
mation having no fundamental physical justification. 
Here we demonstrate that on the contrary KWW is in 
fact a bona fide and respectable relaxation function; ex- 
actly this form of relaxation appears naturally when we 
consider random walks on fractal structures in closed 
spaces of general dimensions. We suggest that the phys- 
ical significance of the relaxation behavior in numerous 
complex systems should be reconsidered in the light of 
this result. 

For random walks on fractal clusters in Euclidean (flat) 
spaces, it is well known that if r(t) is the distance of 
the walker from the starting point after time t, then 
< r 



tion critical exponents |t],|^] whose numerical values are 
known quite accurately in all dimensions |p[-|l4]| . 
For random walks on the surface of a sphere, which is 
a closed surface, the local behavior < r 2 (t) > oc t can 
be shown to lead exactly to an exponential decay of the 
autocorrelation < cos(9(t)) >oc exp(—t/r) where 6(t) is 
the angle between the initial t = position vector of 
the walker and the position vector at time t. With an 
appropriate definition of < cos(6(t)) > this result holds 
for hyperspheres in any dimension. For random walks 
on a percolation fractal inscribed on a closed space with 
the topology of a hypersphere, it was conjectured some 
years ago that the end-to-end autocorrelation function 
should decay as a KWW stretched exponential, with the 
Kohlrausch exponent (3k equal to the flat space percola- 
tion fractal random walk exponent, i.e. (3k = (3rw m am/ 
space dimension [ ^5|Jl6| . This can be understood simply 
in terms of a "fractal" time t RW replacing t in both the 
local < r 2 (t) > and the exponential < cos(8(t)) > expres- 
sions above. This conjecture has been extensively tested 
numerically but only in the extreme case of the very high 
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Here (3rw = d/D 



where d and D are the spectral and fractal dimensions of 
the cluster respectively. For a critical percolation fractal, 
9 = (fi — (3)/v where fi, (3 and v are universal percola- 



dimensional hypercube p6| , p7 22 (the hypercube has the 
same closed space topology as a hypersphere). In the 
present work we have studied numerically the general 
case of random walks on percolation clusters on hyper- 
spherical surfaces for the range of different embedding 
dimensions d running from 3 to 8. 

The simplest case to visualize is the surface of a sphere 
in dimension d = 3. The surface is decorated with small 
disks whose centers are distributed at random. Several 
clusters, made of overlapping disks, can be determined. 
The largest cluster contains more and more disks as the 
total number of disks is increased. Just as there appears 
a percolation cluster (containing a non-infinitesimal frac- 
tion of the total number of disks) for the equivalent sys- 
tem in the two dimensional flat space above a critical 
value of disk concentration, so the largest cluster "per- 
colates" on the surface of the sphere when the number 
of disks is sufficiently large. As an illustration, in fig- 
ure 1, we have represented the two dimensional projec- 
tion of the percolating cluster in the case S/R = 0.01 
(where S is the disk diameter and R the sphere radius). 
Here N p = 179200 disks have been disposed at random; 
only the largest cluster, containing 38130 disks, is repre- 
sented. Note that this cluster spans almost a hemisphere, 
a situation intermediate between a well localized cluster 
(N << N p ) and a cluster spanning uniformly the whole 
surface of the sphere (N >> N p ). 

The position of a given disk center can be defined by the 
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d = 3 coordinates Xi, where the origin is taken as the 
center of the sphere. Imagine now a walker jumping at 
random from one disk center to the center of any disk 
overlapping it. While the values of the coordinates aver- 
aged over many walks < Xi(t) > stay finite for N << N p 
(the walker is localized) they decay exponentially to zero 
for N >> N p , finally losing memory of their initial val- 
ues. Of course, in the limit TV >> N p , < Xi(t) > goes 
to zero since the random walker starts to investigate the 
largest cluster (which fills uniformly the curved space) 
entirely. Our calculations show that right at the perco- 
lation where the largest cluster is fractal, the decay is 
critical, and takes up precisely the stretched exponen- 
tial form, with a (3k exponent equal to (3rw, already 
known for the random- walk on a percolating cluster in a 
two-dimensional flat space. Our calculation is the gener- 
alization of this picture over a wide range of dimensions, 
in particular for dimensions larger than n = 6, where it 
is known that [3rw reaches its mean-field value 1/3. 
Consider the surface of a d-dimensional (hyper)sphere of 
unit radius which can be defined, using Euclidean coor- 
dinates, by Y^i=o x f = 1- This (hyper)surface is a closed 
and curved n-dimensional (n = d— 1) space, S n , on which 
one can define a geodesic distance between two points (1) 
and (2) by 8 — cos^ 1 q where q is the scalar product of the 

end positions, i.e. q — i"i.r 2 = Y^i=o x i x i ■ Of course, 
this distance becomes asymptotically equivalent to the 
Euclidean distance in the limit of distances infinitesi- 
mally small compared to the radius of the (hyper)sphcre 
(which is here set to unity). On this n-dimensional (hy- 
per)surface, identical n-dimensional small (hyper)disks 
(called disks in the following) of diameter d « 1 are dis- 
posed sequentially, the successive disk centers being cho- 
sen at random uniformly on S n . To determine the clus- 
ter structure, as soon as a new disk is added, a search for 
connections with previous disks is performed by checking 
if their center-to-center geodesic distance is smaller than 
S. At each stage, when N disks have been disposed on 
S n , we can define and label the different clusters made 
up of connected disks. In particular one can determine 
the largest cluster containing the greatest number N m 
of disks. Also one can characterize the overall filling of 
the space by a dimensionless parameter 77, ratio between 
the sum of the area of the individual disks and the to- 
tal area of S n Note that this filling parameter is 
larger than the true volume fraction c, because of multi- 
ple counting due to overlaps, and therefore r\ can even- 
tually exceed unity. More precisely it has been shown 
that c = 1 — cxp (—77) ]l8||l9| ]. Percolation occurs for 77 
larger than a threshold value r\ v above which P = N m /N, 
which measures the probability for a given disk to belong 
to the largest cluster, remains non-zero in the thermody- 
namic limit 5 — ► 0. In practice the curve P(jj) exhibits a 
sigmo'idal shape which becomes sharper and sharper at 
i] p as S — > 0. 



Once the largest cluster has been identified, we choose 
one of its constitutive disk centers r(0) at random as a 
starting point. We perform a random walk on the clus- 
ter by performing successive jumps, first from r(0) to 
the center r(l) of any other disk connected to it (chosen 
at random over all its overlapping disks), then iterating 
from r(l) to r(2), etc.. [^TJ After t steps, the "correlation 
function" r(0).r(i) is calculated. This is the scalar prod- 
uct of the two end positions, which is no more than the 
cosine of the geodesic end-to-end distance 8 measured on 
S n . In practice we calculate the quantity q(t) =< cos 8 > 
which has been averaged, for a given number of steps t, 
over N r independent realizations of the largest cluster as 
well as over N s independent choices of the starting point 
on each cluster. Given dimension n and size 5, the be- 
havior of q(t) has been analyzed for different filling values 

7]. 

An example with n — 3 and 6 — 0.02 is shown in figure 
2, where — In \q(t)\ has been plotted as a function of t in 
a log-log plot, after taking an average over independent 
walks starting from N s = 10500 different disks on each 
of the N r = 20 independently generated largest clusters 
(a typical value of the number of disks in the percolating 
cluster is N m ~ 700000). 

If the relaxation function q(t) is strictly a stretched ex- 
ponential, this type of plot produces a straight line of 
slope Pk- For the critical value of 77 (here rj p ~ 0.3452), 
one observes a clear straight line behavior in the numer- 
ical data over (at least) six decades in t (the walk for 
rjp has been extended to illustrate this point). In the 
inset, we show how we estimate rj p and [3 k- The effec- 
tive slope of the log-log plot has been determined by a 
least-square fit within an interval of one-tenth of the total 
range in hit and plotted as a function of Int. The per- 
colation threshold rj p is estimated as the n-value giving 
the widest plateau at large times, and the [3k exponent 
is taken as this plateau value. This procedure has been 
repeated for different values of 6 and for dimensions n 
ranging from 2 to 7 (for n~l it is well known that the 
percolation transition does not occur Jll[). The lowest 
attainable 6 values are mainly determined by the limited 
memory of our computers. In practice the values of N s 
have been chosen to be of the order of two percent of 
N m , and N r has been chosen so as to obtain runs of the 
order of a day on regular high speed PC computers. N s 
varies from few thousands to a few units when decreasing 
5 from about 0.3 to its lowest attainable 6 value. This 
protocol leads to error bars of the order of 0.01 for the 
exponent estimates. Of course the "true" values of the 
percolation thresholds and exponents are obtained by an 
extrapolation to the "thermodynamic" limit 5 — ► 0. 
The numerical results for rjp in 77 = 2 and 3 are in excel- 
lent agreement with accurate estimates from flat space 
calculations [ |l8|]Ts| ], and rj p drops quasi-exponentially 
with n for higher dimensions. The data for [3k are sum- 
marized in figure 3, where [3k has been plotted as a func- 
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tion of 6. On the S — axis of figure 3, we have indicated 
the best estimates of the flat space (3rw with their as- 
sociated error bars, calculated from recent fx, (3 and v 
values available in the literature [PHbj . We note that for 
all n > 6, (3rw = 1/3 exactly, since n = 6 is the upper 
critical dimension for percolation . It is quite remark- 
able that for each dimension n, a simple straight line fit 
of our data goes through the corresponding (3rw value 
to within the numerical error or, at least, extrapolates to 
a value very close to it. 

It should be noticed that the law < cos 8 >oc 
exp — (t/r) l3K not only contains the large time relax- 
ation behavior but also the short time behavior as, af- 
ter expanding both sides for small t and 0, it becomes 
< 8 2 >oc (t/r)P K . Since for short times the random 
walker stays on a n-dimensional surface tangent to the 
(hyper)sphere, one recovers the law < r 2 >cx t@ RW in 
dimension n. This could explain why the stretched ex- 
ponential behavior extends over so large a region of time 
(see figure 2). In practice, as in flat space, at very short 
times there are corrections to scaling due to the discrete 
character of the walk. 

These data can be taken as a clear numerical demon- 
stration that random walks on a fractal cluster inscribed 
on a hypersphere in any dimension lead necessarily to 
a stretched exponential decay of the correlation func- 
tion, with a Kohlrausch exponent (3k equal to the ratio 
d/D of the spectral and fractal dimensions of the cluster. 
Thus the stretched exponential relaxation on a fractal in 
a closed space appears as the precise analogue of the sub- 
linear diffusion on a fractal in a flat space. 
Why is this relevant to relaxation in complex systems ? 
A complex system is made up of many individual ele- 
ments (atoms, molecules, spins, ...etc), all in interaction 
with each other. The total space of all possible configu- 
rations of the whole system is a huge closed space having 
a very high dimension, of the order of the number of ele- 
ments. These spaces are so astronomically large that an 
explicit evaluation of their properties, configuration by 
configuration, is almost impossible except for tiny sys- 
tems. Each configuration has an energy associated to 
it. At finite temperature T only a restricted subset of 
configurations are of low enough energy to be thermo- 
dynamically accessible by the system. In equilibrium at 
T above any ordering temperature, the system is perma- 
nently exploring all this subset of accessible configura- 
tions by successive movements or reorientations of local 
elements of the total system. By definition this equilib- 
rium relaxation can be mapped onto a random walk of 
the point representing the instantaneous configuration of 
the total system within the space of accessible configura- 
tions. Thus the configuration space of a complex system 
can be viewed as a "rough landscape" . In such a scenario, 
the portion of configuration space available to the system 
consists of only a restricted set of tortuous configuration 
space paths. In this situation, when real measurements 



are made, the observed relaxation functions must reflect 
the complex morphology of the available configuration 
space, and so will be slow and non-exponential - typi- 
cally stretched exponential. Now, any relaxation func- 
tion including the KWW function can be represented as 
the sum of an appropriate distribution of elementary ex- 
ponential relaxations. However it is important to note 
that the independent exponentially relaxing elements in 
a complex system are the modes of the whole system not 
the atoms, molecules or spins which are in strong inter- 
action with each other |2^] . It is the morphology of the 
configuration space which determines the mode distribu- 
tion and the form of the relaxation. 
We have just seen that a closed space fractal structure 
necessarily leads to a relaxation process which is exactly 
of stretched exponential type; we suggest that, inversely, 
when a complex system at temperature T is actually ob- 
served to relax with a stretched exponential, it is the 
signature of a fractal morphology of the available config- 
uration space at that temperature. In particular, what- 
ever the microscopic details of the interactions in a glass 
former are, it is generally observed that as the glass tem- 
perature is approached from above, KWW relaxation sets 
in. The implication is that phase space takes up a fractal 
morphology as a consequence of the intrinsic complexity 
associated with glassiness. The image of the glass tran- 
sition which follows is that of a percolation transition 
in phase space [[l5],^2| with the temperature being ana- 
log to the filling parameter. As would be expected from 
the argument given above, the Kohlrausch exponent is 
observed to tend to a limiting value of 1/3 as the glass 
transition is approached in a number of systems, p^-|2r^] 
for instance. 

In summary, we have demonstrated that random walks 
on fractals in closed spaces give stretched exponential 
relaxation, and we suggest that stretched exponential 
relaxation is ubiquitous in nature because configuration 
spaces are fractal in many complex systems. 
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FIG. 1. Two-dimensional projection of a typical perco- 
lating cluster of disks of diameter 8 — 0.01 on the surface 
of a sphere of unit radius. Here N p = 179200 disks have 
been disposed at random; only the largest cluster, containing 
Nm = 38130 disks, is represented. 




FIG. 2. Plot of — In \q(t)\ versus t (averaged over N r = 20 
and N s = 10500) for n = 3, 5 = 0.02 and for different r) 
values: from bottom to top r\ — 0.335, 0.34, 03452, 0.35, 
0.355 (the statistical error bars do not exceed the thickness 
of the lines). In inset the local value of the exponent /3 is 
plotted as a function of In t for the different r/ values in order 
to estimate r\ v and /3k- 
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FIG. 3. Variation of the exponent (3k as a function of the 
(hyper)disk diameter S for dimension n = 2,3,4,5,6 and 7 
(from top to bottom). The solid lines represent least square 
linear fits of the data points. The best estimates of the flat 
space exponent Pbw with their associated error bars [9-14] 
are represented by the stars on the y axis. 
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